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Abstract —  Reference  [1],  which  appears  in  the 
Proceedings  of  the  RADAR’09  Conference  at  Bordeaux  in 
October  2009,  describes  nonlinear  currents  generated  by  an 
electromagnetic  wave  incident  on  a  metal  surface.  This  paper 
supplements  the  results  of  [1]  with  a  full  calculation  of  the 
RCS  per  unit  area  for  second  harmonic  generation.  These 
results  verify  that  the  level  of  such  emission  from  a  metal 
surface  is  anomalously  low. 
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I.  Introduction 

Current  and  future  proposals  to  use  radar  in  situations 
that  extend  beyond  its  traditional  role  as  a  method  of  remote 
sensing,  e.g.,  airport  security  (weapons  detection),  mapping 
(SAR),  and  mine  detection,  are  often  constrained  by  the  need 
to  overcome  certain  limitations  of  its  capabilities.  When  this 
situation  exists,  various  marginal  technologies  need  to  be  re¬ 
examined  periodically  with  a  goal  of  addressing  such  needs. 
One  such  technology  is  harmonic  and  multitone  radar,  in 
which  a  signal  illuminating  a  target  at  one  frequency  creates  a 
response  at  another  frequency.  The  advantage  of  this  type  of 
radar  lies  in  its  frequency  diversity  -  out-of-band  returning 
signals  may  escape  jamming,  clutter  and  multipath  effects  can 
be  mitigated,  etc.  In  addition,  the  receiver  may  be  able  to 
identify  new  properties  of  the  target,  i.e.,  novel  signatures. 

This  paper  will  focus  on  the  last  possibility. 

In  many  respects,  understanding  the  nonlinear 
response  of  a  target,  which  was  the  focus  of  our  paper  Ref.  [1], 
is  primarily  a  problem  of  material  science  rather  than  radar 
technology.  In  Ref.  [1]  we  investigated  the  ability  of  a  signal 
to  induce  nonlinear  currents  at  a  metal  surface,  without 
investigating  the  ability  of  these  currents  to  radiate  energy 
back  into  the  vacuum.  In  this  paper  we  complete  this  task  by 
computing  the  fields  radiated  by  the  second-harmonic 
currents,  and  present  results  for  the  electromagnetic  fields 
radiated  into  vacuum  by  these  currents. 


the  metal  surface  whose  thickness  is  the  metallic  skin 
depth. [2]  Within  the  skin  depth  the  fields  attenuate  rather  than 
propagate,  causing  the  metal  to  reflect  electromagnetic  energy 
efficiently.  At  radar  frequencies  this  high  reflectivity  obeys  the 
Hagen-Rubens  law  [3]: 


9t  =  -l  +  (l-i) 


where  =  107  GHz  is  the  metal  plasma  frequency  and  t  is 

its  conductivity  (Maxwell)  relaxation  time.  In  copper  we  find 

that  co p  =10  GHz  and  C0pT«  320,  so  that  the  reflectivity  at 

10  GHz  is  effectively  -1,  i.e.,  total  reflection  apart  from  a 

correction  of  order  10  E  which  is  in  fact  the  transmittivity 
into  the  metal.  In  addition,  since  the  skin  depth  of  copper  at  10 
GHz  is  only  60  microns,  it  is  clear  that  nonlinear  effects  are 
strongly  localized  at  the  surface. 

In  [1]  we  discussed  a  well-established  model  for  the 
metal  electron  “gas”  based  on  plasma  physics  [4],  where  the 
motion  of  mobile  electrons  obeys  the  equations  of  motion  for 
a  charged  fluid  (i.e.,  a  plasma).  Because  the  radar  signal  is 
never  very  strong  at  a  target,  it  is  legitimate  to  solve  these 
equations  in  a  perturbation  series.  The  various  terms  in  the 
perturbation  series  generate  new  frequencies  in  the  induced 
target  current,  which  radiate  energy  back  to  the  receiver. 


III.  ILLUMINATION  IN  THE  LINEAR  REGIME 

Let  us  begin  with  the  case  of  a  weak  EM  field 
incident  on  a  metal  surface.  .  It  is  well  known  [5]  that  in  order 
to  specify  the  problem  of  off-normal  incidence  of  an  EM  wave 
on  a  surface,  we  must  first  define  a  plane  of  incidence,  which 
contains  both  the  wave  vector  of  the  incident  wave  and  the 
surface  normal.  Then  an  arbitrary  wave  polarization  will  be  a 
linear  combination  of  two  components: 


II.  Bare  metal  surfaces 

As  discussed  in  [1],  electromagnetic  fields  are 
difficult  to  create  within  a  metal  due  to  screening  by  its  mobile 
electrons,  which  confines  the  incident  fields  to  a  region  near 


(1)  Magnetic  field  in  the  plane  of  incidence  and  electric 
field  perpendicular  to  it  (s  -  polarization) 

(2)  Electric  field  in  the  plane  of  incidence  and  magnetic 
field  perpendicular  to  it  (/?  —  polarization) 
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The  geometry  of  this  problem  is  shown  in  Fig.  1  for  a  p  — 
polarized  wave.  This  polarization  tends  to  have  the  strongest 
interaction  with  the  metal  since  the  electric  field  generally  has 
a  component  perpendicular  to  the  surface.  For  this  case  the 
linear  EM  fields  outside  the  metal  (z  >  0)  are 


E=  Er. 


-ikz  cos0g_  +3{eikz  cos0^ 


B  = 


]Eo(e~‘ 


- ikz  cos  0  _  <ftjkz  cos  0 


/Axsin0  —/cot 
e  e 


/Ax  sin  0  -/cot 
e  e 


z 


Figure  1 :  p  -  polarized  incident  wave  geometry;  the  unit 
vectors  e;,  are  polarization  (electric-field)  vectors.  Vectors 

are  written  in  rectangular  component  form,  i.e.,  (Ax  ,  Az) 

where  Eq  is  the  incident  electric  field,  is  the  reflection 

coefficient  of  the  surface,  k  =  col  c ,  and  0  is  the  angle  of 
incidence  measured  form  the  surface  normal.  Within  the 
metal,  the  wave  properties  are  modified  by  the  presence  of 
mobile  charges.  As  discussed  in  [1],  combining  the  plasma 
equations  with  the  Maxwell  equations  leads  to  both  electrical 
and  mechanical  degrees  of  freedom  in  the  metal,  which  are 
coupled  in  such  a  way  that  the  electrical  oscillations  in  the 
metal  fall  into  two  categories  [6] : 

(1)  Electromagnetic  oscillations  with  nonzero  magnetic 
field  and  zero  density  oscillations  (referred  to  as 
“transverse  plasmons”),  and 

(2)  Electrostatic  oscillations  with  zero  magnetic  field  and 

nonzero  density  oscillations  (referred  to  as 
“longitudinal  plasmons”) 


Maxwell  boundary  conditions,  and  leads  to  dispersive  surface 
waves  (called  “surface  plasmons”)  that  have  been  studied 
extensively. 

Within  the  metal  (z  <  0)  the  fields  take  the  form 


E  = 


B  = 


Ae 


■  P?  ( *  /A  sin  0  A  -Oz(ik  sin0 

V  -I - -7  -4-  Afz?  - 


icP 


x  H - z  +  Be 

P 


Pz  - 


/Ax sin 0  -/co/ 
e  e 


CO  CO+- 
V  x 


f  Ae  y 


g/AxsinO^-Zco/ 


where  A  and  B  are  determined  by  the  continuity  conditions 
plus  the  extra  boundary  condition.  The  angle-dependent 
attenuation  exponents  are  given  by 


1 

CO 

2 

i 

c 

COH - 

X 

(»p  -  co 


(1) 


Q2  =  k2  sin2  9  +  — y|  ®p2 


where 


(Op  = 


n0e 

e  Q>n 


co+ 


fi  (-  2  \1/3 

V  ~  -l3*  "0  j 


are  the  “plasma  frequency”  and  Fermi  velocity  of  the  metal. [8] 
Some  manipulation  of  the  fluid  equations  shows  that  the 
current  is  given  by 


—/(OH — 

X 

Using  the  expression  for  E  inside  the  metal,  we  find  that 
-pz 

jyz  =  kriC.Be  "  ,  so  that  the  boundary  condition 

/|  ^  (z  =  0)  =  0=>A?  =  0.  The  constant  A  can  then  be  written 

3fiQ ,  thereby  defining  a  transmission  coefficient  3  along 

with  dt.  Applying  the  electromagnetic  boundary  conditions 
then  gives  us 


co  p1 2E - V  (V  •  eJ 

VF 


Dissipation  causes  such  oscillations  to  be  attenuating  in  time 
when  free  and  evanescent  in  space  when  driven  by  a  field 
from  outside  the  metal.  The  additional  mechanical  degrees  of 
freedom  (electron  density  and  velocity)  imply  that  the  usual 
Maxwell  boundary  conditions,  i.e.,  continuity  of  tangential  E 
and  B  ,  are  insufficient  to  define  the  problem  in  the  presence 
of  the  surface,  and  require  an  additional  boundary  condition, 
which  in  the  metal-physics  literature  is  taken  to  be  the 
condition  that  there  be  no  electron  velocity  perpendicular  to 
the  surface.  [7]  This  condition  is  independent  of  the  usual 
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These  expressions  solve  the  linear  problem  completely,  which 
initializes  the  perturbation  calculation. 

IV.  Nonlinear  currents 


The  second-order  equation  for  the  electric  field  is  an 
inhomogeneous  version  of  the  linear  equation  of  the  form 


(A+rl 
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where  S  =  5^  +  are  source  terms  arising  from 

the  four  nonlinear  processes  discussed  in  [1].  For  an  incident 
single-frequency  p  -  polarized  plane  wave  we  can  introduce 

the  usual  complex  notation  in  which  E  =  Re  E  and  S  =  Re  E , 
in  which  case  these  inhomogeneous  currents  become: 
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Since  the  x  -  dependence  of  the  fields  must  match  that  of  the 
currents,  the  second  order  equation  reduces  to  two  coupled 
ordinary  linear  differential  equations  for  the  electric  field 
components  E 0  and  E^  ^  : 
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where  K  =  k  sin  0  and  the  inhomogeneous  terms  are  known 
and  rather  simple  functions  ofz.  Since  they  are  quadratic 
combinations  of  the  linear  quantities,  which  in  turn  are 
decaying  exponentials  into  the  metal  bulk,  we  can  write 


/Cl)£' 

m 


,r,  2n2  2iKx  —2iS)t 

where  the  coefficients  in  the  large  brackets  are  given  in  the 
Appendix.  This  in  turn  implies  that  the  solution  to  the  field 
equations  will  have  a  particular  part  and  a  homogeneous  part, 
where  the  former  must  have  the  form 


A  = 


2ro 

„  i 
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3c 
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2  2 
c  c 


h  =  AP 


The  matrix  M(n)  arises  from  the  linear  (vector)  operator  on 

the  left  side  of  Eqs.  (2a)  and  (2b).  Such  operators  always 
appear  in  this  type  of  perturbation  expansion. 


For  the  homogeneous  solution  we  must  solve  another 
boundary  value  problem.  Equations  (2a)  and  (2b)  describe  the 
electric  field  inside  the  metal;  outside  in  the  vacuum  this  field 
has  the  form 
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EI ere  the  decay  terms  P, L,  and  are  given  by  Eq.  (1)  with  the 
replacements  CO  — >  2co  and  K  — >  2 K  .  Note  that  P '  ^  2 P  and 
Q 1  jt  2 Q  in  general.  Matching  the  exponentials  gives  a  set  of 
matrix  equations  for  the  particular  solution  fields  of  the  form 


where  C  is  yet  another  unknown.  Clearly  it  is  the  most 
important  quantity  to  determine,  since  it  contains  information 
about  the  fields  leaving  the  metal  and  headed  for  the  receive 
antenna.  Note  that  it  is  an  outgoing  wave  in  this  geometry  with 
no  incoming  wave,  which  is  what  we  expect  for  a  field 
radiation  by  the  nonlinear  currents.  Matching  tangential 
electric  and  magnetic  fields,  and  again  ensuring  that  the 
normal  current  vanishes  at  the  surface,  we  derive  expressions 
for  all  the  unknowns.  The  most  important  one,  C  ,  is  given  by 
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It  is  worth  noting  that  the  denominator  A  has  a  complex  pole 
at  a  frequency  where  =  -ip 2  ’  '•c  -  vvbcrc 


emission  of  power  at  the  target  surface  is  more  than  320  dB 
down  from  the  incident  power,  in  Fig.  2  we  plot  the  magnitude 
of  C  for  copper  versus  angle  of  incidence  of  the  fundamental. 
Careful  algebra  shows  that  C  =  0  for  normal  emission. 


increasing  to  ~  10  0  at  nearly  grazing  emission.  Following 


Ishimaru’s  definition  [10],  the  factor 


C3^ 


is  the 


(power-dependent)  RCS  per  unit  area  for  the  second-harmonic 
signal. 


A  third  factor  that  inhibits  the  use  of  nonlinear  radar 
in  general  arises  from  modification  of  the  radar  equation  for 
nonlinear  radar,  which  we  discussed  in  detail  in  Ref.  [1 1]. 
According  to  the  usual  arguments  [12],  a  field  radiated  by  an 
transmitting  antenna  and  incident  on  a  target  a  distance  R  from 

the  latter  suffers  a  loss  proportional  to  R  '  (which  in  our  case 
here  is  contained  in  the  factor  ).  When  the  target  is 

nonlinear,  this  field  must  be  squared  in  the  second  harmonic 
_2 

current,  producing  an  R  dependence  in  the  radiated  field 

/?  .  Counting  the  return-trip  loss  of  another  factor  of 

2vac 

R  1 ,  it  is  clear  that  the  field  received  by  the  antenna 
(monostatic  case)  is  down  by  R  ,  so  the  power  is  down  by 
A-6  [11]. 


(P2+p2E,)e2-«2(l-E,)  =  0  ■  (3) 

This  expression  was  derived  by  J.  Crowell  and  R.  Ritchie  in 
1968  [9]  as  the  dispersion  relation  for  surface  plasmons  at  a 
metal  surface.  The  existence  of  such  poles  in  the  linear 
reflectivity  of  a  dielectric  and  their  association  with  surface 
waves  is  well  known;  the  appearance  of  this  quantity  Eq.  (3) 
indicates  that  whenever  the  second  harmonic  can  give  rise  to 
this  type  of  surface  plasmon  the  value  of  C  will  be  large. 
However,  since  the  frequencies  where  this  can  occur  are  in  the 
optical  range,  this  result  is  of  little  interest  in  the  present 
context 

V.  Properties  of  second  harmonic  emission 

The  weakness  of  the  second-harmonic  emission  from 
the  metal  surface  can  be  traced  to  two  factors.  The  first  is  the 

2 

presence  of  the  squared  transmission  coefficient  3  in  the 
field  amplitude,  which  implies  a  power  emission  proportional 

to  |3|  .  Because  |3|  in  copper  is  only  ~  10  ^  ,  we  have  an 

immediate  loss  of  160  dB  in  the  amount  of  incident  power 
converted  to  second-harmonic  power.  The  second  factor  is  the 
surprising  inefficiency  of  the  nonlinear  mechanisms  in  the 
metal,  as  measured  by  the  extreme  smallness  of  the  coefficient 
C  ,  which  when  evaluated  numerically  for  copper  gives  rise  to 

—8 

an  additional  decrease  of  ~  10  in  the  amplitude,  so  the  net 


lOLogjQ  |C|  ,  dB 


FIG.  2.  Second-harmonic  RCS 


It  is  easily  shown  that  the  second  harmonic  power  is 
radiated  collinearly  with  the  linearly  reflected  signal.  This 
follows  from  the  replacements  CG  — »  2co  and  K  — >  2 K  in  the 
various  expressions  given  above,  which  implies  that  the  angle 
of  emission  is  identical  with  the  incident  angle.  This  is 
ultimately  due  to  the  evanescent  nature  of  the  fields  inside  the 
metal. 

VI.  Conclusion 

We  have  shown  that  the  strong  suppression  of  second- 
harmonic  emission  from  a  metal  surface  is  a  consequence  of 
two  things:  the  inability  of  the  external  fields  to  penetrate  the 
metal  surface  due  to  linear  reflectivity,  and  the  inefficiency  of 


978-1  -4244-5813-4/ 1 0/S26.00  ©2010  IEEE 


000389 


nonlinear  processes  in  the  metal  associated  the  with  free- 
carrier  (i.e.,  plasma)  mechanisms.  These  results  exacerbate 
difficulties  we  have  discussed  in  Ref.  [11],  and  suggest  that 
the  uses  of  nonlinear  radar  may  be  few.  However,  it  is  worth 
noting  that  because  metals  are  highly  detectable  to  linear 
radar,  there  are  ongoing  efforts  to  render  objects  invisible  by 
making  them  nonmetallic,  e.g.,  plastic  guns.  Such  objects  may 
be  much  more  visible  to  probes  based  on  the  nonlinear 
techniques. 


[11]  F.  Crowne,  C.  Fazi,  and  M.  Ressler,  Proc.  Triservice  Radar  Conf., 
Boulder  Colorado,  June  22-26  2009. 

[12]  S.  A.  Hovanessian,  Radar  System  Design  and  Analysis  (Artech  House 
Inc.,  Dedham  MA,  1984),  Ch.  1. 


VII.  Appendix 

The  coefficients  of  the  inhomogeneous  terms  that  appear 
in  Eqs.  (4a)  and  (4b)  are: 


N  =  1  ^ 

\H z 


c2P  c 2 
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where 


A  =  -^  .  E  =  1  — 


<D+- 


1  (0+- 


The  remaining  quantities  are  defined  in  the  main  text. 
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